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Abstract. A point particle of mass fi moving on a geodesic creates a perturbation 
hab, of the spacetime metric gab, that diverges at the particle. Simple expressions are 
given for the singular fj./r part of hab and its tidal distortion caused by the spacetime. 
This singular part /i^j, is described in different coordinate systems and in different 
CN ' gauges. Subtracting /i^j, from hab leaves a regular remainder h^f,. The self-force on the 

^ I particle from its own gravitational field adjusts the world line at 0(/i) to be a geodesic 

of Qab + h^b'j this adjustment includes all of the effects of radiation reaction. For the 
. case that the particle is a small non-rotating black hole, we give a uniformly valid 

approximation to a solution of the Einstein equations, with a remainder of 0(/i^) as 
/i^ 0. 



, An example presents the actual steps involved in a self-force calculation. Gauge 

' freedom introduces ambiguity in perturbation analysis. However, physically interesting 

problems avoid this ambiguity. 
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■ 1. Introduction 



A description of motion always entails approximations and abstractions. The motion 
of a small black hole through spacetime is clearly not a geodesic of the actual, 
physical spacetime geometry. After all, the "center" of a black hole is inside the event 
horizon, where the geometry is unknown. Nevertheless, if the mass of the hole is 
sufficiently small in comparison with a length scale of spacetime, then the motion is 
approximately geodesic on an abstract spacetime which is described as "spacetime with 
the gravitational field of the black hole removed" . Much of this manuscript focuses upon 
the meaning of this last phrase. 

In general relativity, an object of infinitesimal mass and size moves through a 
background spacetime along a geodesic. If the particle has a small but finite mass fi then 
its world line F deviates from a geodesic of the background by an amount proportional 
to fi. This deviation is sometimes described as resulting from the "self-force" of the 
particle's own gravitational field acting upon itself and includes the effects which are 
often referred to as radiation reaction. 
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In the literature the phrase "gravitational self-force" often refers to precisely the 
right hand side of given below. As emphasized by Barack and Ori 1], the value 
of this quantity depends upon the gauge being used (see section 19 .HI) and is, thus, 
ambiguous. In this manuscript the phrase "gravitational self-force" is only used in an 
imprecise, generic way to describe any of the effects upon an object's motion which are 
proportional to its own mass. 



1.1. Newtonian self-force example 

Newtonian gravity presents an elementary example of a self-force effect |2j . A small mass 
in a circular orbit of radius R about a more massive companion M has an angular 
frequency Q given by 

= M. (1) 

R^l+fi/M)^- ^' 

When fi is infinitesimal, the large mass M does not move, the radius of the orbit R 
is equal to the separation between the masses and = M/R^. However when n is 
finite but still small, both masses orbit their common center of mass with a separation 
of + /i/M), and the angular frequency is as given in (^. The finite /i influences the 
motion of M, which influences the gravitational field within which n moves. This back 
action of fi upon its own motion is the hallmark of a self-force, and the n dependence 
of (0) is properly described as a Newtonian self- force effect. When n is much less than 
M, an expansion of provides 

2/i/M + 0(/xVM')] . (2) 

The finite mass ratio fi/M changes the orbital frequency by a fractional amount 
An _ fi 

In this manuscript we describe any such 0{fi/M) effect on the motion as being a 
"gravitational self- force" effect. Below, we see that the self- force effects for gravity 
include all of the consequences of what is often referred to as "radiation reaction." 
However, we also see that a local observer, near deep inside the wave-zone and not 
privy to global spacetime information, is unable to distinguish radiation reaction and 
the other parts of the gravitational self-force from pure geodesic motion, at this level of 
approximation. 



(3) 



1.2. Electromagnetic radiation reaction in flat spacetime 

The Lorentz force law 

F = g(E + vxB) (4) 

describes the interaction of a point charge q with an electromagnetic field. In an 
elementary electricity and magnetism course, it is implicit that g's own electromagnetic 
field is not to be included on the right hand side — after all for a point charge E is infinite 
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at the very location where it is to be evaluated in (0]). Thus, the electromagnetic field of 
(jU is an "external" field, whose source might be, say, the parallel plates of a capacitor 
but does not include the charge q itself. 

Abraham and Lorentz first derived the radiation reaction force on a point charge 



in terms of the changing acceleration of q. This equation may be interpreted in a 
perturbative sense: Let q have a small mass and be oscillating on the end of a spring. 
At lowest order in the perturbation, q executes simple harmonic motion. At first order 
in the perturbation, the right hand side of (0) is evaluated by use of a v consistent 
with the harmonic motion. The resulting F is a small damping force which removes 
energy from the system at just the proper rate to account for the outward energy fiux 
of radiation. 

A great value of (0) resides in its elementary use by a theorist to calculate the 
radiation reaction force. 

A drawback of © is the apparent obscuration of the root cause of this force. 
Charges interact with electromagnetic fields via Yet, no electromagnetic field is 
present in Imagine a local observer extremely close to q, deep within the wave zone, 
and with a length scale very much smaller than that associated with the oscillations. 
This observer correctly interprets the majority of the acceleration of q as resulting from 
the coupling to the spring. The local observer is unaware of the radiation — a non- 
local concept; yet, he must explain the deviation from pure harmonic motion resulting 
from Frad as a consequence of the interaction of q with some external field via (j3}. 
The Abraham-Lorentz analysis correctly calculates the electromagnetic self-force. But 
it does not explain this self-force in terms of the charge interacting with an external 
electromagnetic field. 

Dirac jl] removes this drawback by providing an interpretation of (jSj) as a direct 
consequence of dH), with the electromagnetic field on the right hand side being an 
external field of indeterminate origin to the local observer. Dirac uses the conservation 
of the electromagnetic stress-energy tensor in a world-tube surrounding q, and ultimately 
takes the limit of vanishing radius of the world-tube. One consequence of his analysis 
is that the half-advanced plus half- retarded field Ff^, = l{F^f + F^"') of q exerts no 
force on q itself, even though the field is formally singular in the point charge limit. We 
call the actual field -F^fe*, and the remainder = F^^* — F^^^ is a vacuum solution of 
Maxwell's equations. F^ substituted into the right hand side of (jl)) yields ©, as shown 
by Dirac. 

A local observer measures the electromagnetic field in the vicinity of q, but with 
no information regarding boundary conditions or distant radiation, he can make no 
conclusions as to the detailed cause or source of the field. However, in the perturbative 
sense described above, the observer can calculate the singular field Fj, in the vicinity of 



[HI 




(5) 
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q. He can subtract this singular field F^^, from the actual, measured field F^^^ to obtain 

F^, = FX - Fl. (6) 

The charge q then interacts with the resulting regular source-free electromagnetic field 
F^ via with a resulting small perturbation in its motion. Thus, a local observer 
naturally explains the damping of the harmonic motion as a consequence of q interacting 
with an external, locally source-free field F^. However, with no global information 
regarding boundary conditions he would not be able to determine the source or cause 
of this external field. In particular the local observer would see no phenomenon which 
he would be compelled to describe as radiation reaction. 



1.3. Electromagnetic radiation reaction in curved spacetime 

DeWitt and Brehme's pioneering analysis of electromagnetic radiation reaction in 
curved spacetime follows Dirac's approach and also uses the conservation of energy in a 
world-tube to determine the force on a point charge. Their results reduce to Dirac's in 
the fiat spacetime limit. However, DeWitt and Brehme find that generally ^{^ab^ + F^b^) 
does, in fact, exert a force on the charge in curved spacetime. After its removal from 
the actual field, the remainder does not serve as the electromagnetic field on the right 
hand side of (j3)) for calculating a radiation reaction force. 

To simplify the remainder of this introduction we, henceforth, assume that the 
charge is in free fall in curved spacetime — the charge would move along a geodesic 
except for interaction with its own electromagnetic field; there are no springs attached. 

DeWitt and Brehme use the Lorenz gauge, Va^° = 0, and a Hadamard expansion 
to break the Green's function into the "direct" and "tail" parts with the vector potential 

^ + Af^ (7) 

The direct part of the retarded Green's function has support only on the past null 
cone, and the tail part has support only inside the past null cone. They find that the 
electromagnetic self-force can be described as a consequence of the particle interacting 
just with Al^"^, 

= g^7"'(Vc4"' - VbAfV- (8) 

This expression, like has the great value that it can be used to calculate an 
electromagnetic self-force, but it shares the drawback that it does not explain the self- 
force in terms of a locally measurable, source free solution of the Maxwell equations. In 
fact A^^'^ is not in any sense a solution of the electromagnetic field equation 

The details of the Hadamard expansion reveal that if A^^^jj were inserted into the left hand 
side here, it would yield a phantom J^" throughout a neighborhood of the charge. There 
would be no other evidence for the existence of this J^^^^. Further, if {Rab — \gabR)u^ 7^ 0, 
then A**^'' is not differentiable at the particle and some version of averaging around the 
charge is required to compute the self-force. A'^^^-y is a valuable mathematical construct 
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which may be used to calculate the self- force from (jSI), but it is not associated with 
an actual electromagnetic field. We conclude that the DeWitt-Brehme construction 
correctly calculates the electromagnetic self-force. But it does not explain the self-force 
in terms of the charge interacting with an external electromagnetic field. 

A modification [Hj of the DeWitt and Brehme analysis has rectified this shortcoming. 
The actual vector potential may be decomposed as 

^Al + A^, (10) 

where and A^ are, in fact, solutions of Maxwell's equations in a neighborhood of 
q: A^ has only the charge q as its source, while A^ is a vacuum solution. Further, (jH]) 
yields the same force whether A^ or A^^^^ is inserted on the right hand side, after the 
possible lack of differentiability of Ajf'' is handled properly. 

One nuance of the decomposition into S- and R-fields, is that the Green's function 
for the S-field has support at the advanced and retarded times, just as in the flat- 
spacetime example, above. But it also has support at the events between the retarded 
and advanced times — these have a spacelike separation with the field point. 

The "S" and "R" decomposition provides a local observer in curved spacetime with 
the ability to measure the actual electromagnetic field F^^^^ in a neighborhood of q. He 
can make no conclusions as to the detailed cause or source of the field. However, in the 
perturbative sense described above, the observer can calculate in a neighborhood 
of q based upon its approximate geodesic motion. He can then subtract this singular 
field F^fj from the actual, measured field -F^^*. The charge q then interacts with the 
resulting regular source-free electromagnetic field F^ via (0} or (jH)) with a resulting 
small perturbation of its geodesic motion. Thus, a local observer naturally explains 
the lack of geodesic motion of a charge g as a consequence of q interacting with an 
external, locally source-free electromagnetic field. However, with no global information 
regarding boundary conditions he is not able to determine the source or cause of this 
external field. In particular, at this level of approximation the local observer sees no 
phenomenon which he would be compelled to describe as radiation reaction. 

Gravitational self-force 

The treatment of gravitational radiation reaction and self-force, in terms of Green's 
functions, are formally very similar to that just described for the electromagnetic field. 

In some circumstances the gravitational field may be considered to have an effective 
stress-energy tensor consisting of terms which are quadratic in the derivatives of the 
metric. Mino, Sasaki and Tanaka [Zj follow the DeWitt-Brehme approach, but with 
this gravitational stress energy tensor. Ultimately, they conclude that the motion of a 
point mass fi satisfies 

fiu'Vtu^ = -fxig""' + u^u'Xu'^iVXf - ^V,/i*f ). (11) 

In an independent analysis within the same paper, they treat as a small black hole 
moving in an external universe and use a general matched asymptotic expansion to 
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arrive at the same conclusion. In this latter approach, the metric of the black hole 
is considered to be perturbed by the external universe through which it is moving. 
Simultaneously, the metric of the external universe is considered to be perturbed by the 
small mass fi moving through it. Others have used matched asymptotic expansions to 
describe the motion of a small black in an external universe jHl El UHl UH 1121 HHj, but 
the connection between such results and radiation reaction appears not to have been 
made before reference [7j. 

Quinn and Wald [Hj use an axiom based analysis of the gravitational self-force and 
also arrive at 

The form of equation (jllll is equivalent, through first order in to the geodesic 
equation for the metric Qab + h^iX^. From one perspective then is the gravitational 
equivalent of (jHl). Equation (fTT| . like has the great value that it can be used to 
calculate a gravitational self-force, but it shares the drawback that it does not explain 
the gravitational self-force in terms of geodesic motion in a locally measurable, source 
free solution of the Einstein equations. In fact, /i*'^'^ is not in any sense a solution of the 
perturbed Einstein equation, given below in (fT3|) . 

The details of the Hadamard expansion reveal that if h^^j^^ were inserted into the 
left hand side of (fT^ . it would yield a phantom stress-energy tensor T^j^'^ throughout 
a neighborhood of fi. There would be no other evidence for the existence of this T*^''. 
Further, when RacbdU^W^ 7^ 0, then /i*^'' is not even differentiable at the particle; although 
details reveal that averaging around the particle is not required to compute the self-force 
with (|TT|l . h}^^^ is a valuable mathematical construct which may be used to calculate 
the self-force from ffTTj) . but it is not associated with an actual gravitational field. We 
conclude that the Mino, Sasaki and Tanaka and the Quinn and Wald constructions 
correctly calculate the gravitational self-force. But they do not explain the self-force in 
terms of geodesic motion in an external gravitational field. 

A modification jH] of the analysis involving /i*^'^ has rectified this shortcoming. The 
actual metric perturbation may be decomposed as 



where h^-^ and h^^ are, in fact, solutions of the perturbed Einstein equations (fT!?|l in a 
neighborhood of /x: h^^ has only the mass /x as its source, while h^^^ is a vacuum solution. 
Further, (|TT|) yields the same force whether /i^^ or h}"^^^ is inserted on the right hand side. 

Earlier ^2j; asymptotic matching was used to find an explicit expression for the 
leading terms in an expansion of /i^^ in powers of the distance away from /i. Further, 
it was also shown that h^^^ = h^f — /i^^ was at least C^, with the given terms of the 
expansion for /i^^, and that n necessarily followed a geodesic of gab + h^b ^'^ terms 
of 0(/i^/7^^), where 71 is a length scale of the background geometry. However, at that 
time it was erroneously claimed that the /i^^ field was identical to /i*^'' because both 
led to the same equation of motion — namely geodesic motion in gab + h^h- during 
a failing effort to demonstrate directly this equivalence that the important differences 
between the pair h^^ and /t^^ and the pair h'^^,^ and h^^f^^ as possible solutions of the 




+ h, 



R 

abi 



(12) 
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perturbed Einstein equations were discovered [Uj. 

A small mass /i moves through a background geometry gab along a world line F. At 
the lowest order in a perturbative sense, F is a geodesic. The Newtonian example given 
in section [TTT] implies that F deviates from geodesic motion in gat by O^ji/lZ) — it is this 
deviation in which we are interested. 

A local observer in curved spacetime has the ability to measure the actual metric 
g^ in a neighborhood of /i. In a perturbative sense, the observer can calculate /i^^ 
in a neighborhood of /i based upon its approximately geodesic motion. He can then 
subtract this singular field /i^j, from the actual, measured field g^. The mass will 
be observed to move along a geodesic of g''^ — h^^^ = gab + h^b- Thus, a local observer 
sees geodesic motion of /i in the metric gab + h^b^ which is a vacuum solution of the 
Einstein equations, up to a remainder of 0(/x^) in a neighborhood of fi. With no global 
information regarding, say, the original background metric gab, he would be unable to 
make any measurement which would distinguish the separate parts gab and /i^j, which 
together make up the metric through which /z is moving on a geodesic. At this level of 
approximation the local observer sees only geodesic motion and no phenomenon which 
he would be compelled to describe as radiation reaction. 

1.5. Outline 

Perturbation analysis, described in|2l is the heart of the self-force formalism. A variety of 
locally inertial coordinate systems are identified in|21 Some of the ensuing mathematics 
is simplified by use of notation, introduced in|3J which is convenient for describing vector 
and tensor harmonics in a spherically symmetric geometry. 

Sections|31|7|describe the metric in the neighborhood of a small black hole as it moves 
through spacetime and provide an identification of the singular "S-part" of a particle's 
gravitational field, which exerts no force on the particle, itself. The remaining "R-part" 
of the particle's gravitational field is then seen to be responsible for the gravitational 
self-force in |S] The confusion caused by the gauge freedom inherent in the R-part is 
summarized in El 

An example of a point mass in a circular orbit about a Schwarzschild black hole 
reveals, in section[TUl how the difficulty of gauge dependence may be handled in carefully 
defined circumstances. Future prospects for gravitational self-force calculations are 
discussed inlTTl 

1.6. Conventions and notation 

Conventions and notation are described here and again in context below. The indices 
a, b, c. . . are spacetime indices lowered and raised with the metric gab and its inverse; 
the derivative operator compatible with gab is Va. The metric of fiat Minkowskii space 
is rjab- The indices i, j, k, I, p, q are always used as spatial indices and are raised and 
lowered with the fiat three-metric fij. hi is a unit radial vector in fiat space. 
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Indices A, B, . . . are used to denote vector or tensor components which are tangent 
to a two-sphere in spherically symmetric geometries, especially those which are generated 
by "potential" functions as described in section HI Spatial, symmetric trace-free tensors 
such as £ij or Bijk represent the external gravitational multipole moments, when the 
gravitational field field is expanded in a locally inertial coordinate system. The symbols 

5 and B always refer to the even and odd parity moments, respectively. The scalars 
£^'^^ = SijU^W and B^^"* = Bij^n'^Wn'' , for examples, represent linear combinations of the 
i = 2 and i = 3 spherical harmonics, respectively, which depend only upon the angles 

6 and in the usual Schwarzschild coordinates, and are independent of t and r. the 
superscript (2) denotes the value of i. 

A small particle of mass fi moves along a world line F parameterized by the proper 
time s. p is an event on T. 7^ is a representative length scale associated with a geodesic F 
of spacetime — 71 is the smallest of the radius of curvature, the scale of inhomogeneities, 
and the time scale for changes in curvature along F. We use h^^^o represent the singular 
source field, while is an approximation to h^f, based upon an asymptotic expansion. 



2. First order perturbation analysis 

Perturbation analysis provides the framework for an understanding of the self-force and 
radiation reaction on an object of small mass and size in general relativity. This begins 
with a background spacetime metric Qab which is a vacuum solution of the Einstein 
equations Gab{g) = 0. An object of small mass then disturbs the geometry by an 
amount hat = 0(/i) which is governed by the perturbed Einstein equations with the 
stress-energy tensor Tab = 0(/i) of the object being the source, 

Eabih) = -8nTat + 0ifi^). (13) 

Here Eab{h) is the linear, second order differential operator on symmetric, two- indexed 
tensors schematically defined by 

Ea,{h) ^ -^/led, (14) 

and Gab is the Einstein tensor of gab-, so that 

2Eab{h) = V^hab + Va^bh " 2V^aV'h), 

+ 2Ra\%d + gabiV'VKd - V^/l), (15) 

with h = habg""^ and Va and Ra'^b'^ being the derivative operator and Riemann tensor of 
gab- If hab is a solution of (|T!^ then it follows from (fT^ that gab + hab is an approximate 
solution of the Einstein equations with source T^^, 

Gab{g + h)=MTab + 0{lx'). (16) 

The Bianchi identity implies that 

V^'Eabih) = (17) 
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for any symmetric tensor hat', this is discussed in Appendix A Thus, an integrabihty 
condition for (fT^ is that the stress-energy tensor Tab be conserved in the background 
geometry gab, 

VXfe = 0(/i2). (18) 

Perturbation analysis at the second order is no more difficuh formally than at 
the first. But the integrabihty condition for the second order equations is that Tab be 
conserved not in the background geometry, but in the first order perturbed geometry. 
Thus, before solving the second order equations, it is necessary to change the stress- 
energy tensor in a way which is dependent upon the first order metric perturbations. 
This modification to Tab is said to result from the "self-force" on the object from its 
own gravitational field and includes the dissipative effects of what is often referred to 
as "radiation reaction" as well as other nonlinear aspects of general relativity. This 
modification to Tab is 0(/i^) because Tab itself is 0(/i). 

A description of general, nth order perturbation analysis is given in Appendix B 



The procedure is similar to that just outlined. The stress-energy tensor must be 
conserved with the metric fi'^b"^'' order to solve the nth order perturbed Einstein 
equation ()B.4|) for h^ff^ . In an implementation, the task then alternates between solving 
the equations of motion for the stress-energy tensor and solving the perturbed Einstein 
equation for the metric perturbation. Similar alternation of focus between the equations 
of motion and the field equations is present in post-Newtonian analyses. 

For many interesting situations the object is much smaller than the length scale of 
the geometry through which it moves. We expect, then, that the detailed structure of 
the source should be unimportant in determining its subsequent motion. 

To focus on those details of the self-force which are independent of the object's 
structure we first attempt to model the object by an abstract point particle with no spin 
angular momentum or internal structure. The stress-energy tensor of a point particle is 

/oo n.a b 
^^(5^(x'^-X'^(s))ds (19) 
■oo V-9 

where X°'{s) describes the world line F of the particle in some coordinate system as a 
function of the proper time s along the world line. 

The naive replacement of a small object by a delta-function distribution for the 
stress-energy tensor is satisfactory at first order in the perturbation analysis. The 
integrabihty condition ()18|) requires the conservation of the perturbing stress-energy 
tensor. For a point particle this implies that the world line F of the particle is 
an approximate geodesic of the background metric gab, with u'^V aU^ = 0(/i) (c/ 
[Appendix CD . The solution of ()13|1 is formally straightforward, even for a distribution 
valued source. This procedure has been used many times to study the emission of 
gravitational waves by a point mass orbiting a black hole I16| ITTIj. 

A difficulty appears with the second order integrabihty condition ()B.10|) . with n = 2. 
This condition seems to require that the particle move along a geodesic of gab + hab- But 
hab is singular precisely at the location of the particle. To rectify this situation we look 
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for a method to identify and to remove the singular part h^^ of the point particle's metric 
perturbation and, thus, to find the remaining /i^^. We would then have the expectation 
that the point particle would move along a geodesic of the abstract, perturbed geometry 
9ab + h^b- 

To avoid the singularity in hab, we replace the point particle abstraction by a small 
Schwarzschild black hole. The difficulty caused by the formal singularity is replaced by 
the requirement of boundary conditions at the event horizon. Following Mino, Sasaki 
and Tanaka [Zj , in section IHl we use a matched asymptotic expansion to demonstrate 
how the 0{fi) self-force adjusts the world line of the particle. For a small black hole 
moving in an external spacetime, the solution of the Einstein equations divides into two 
overlapping parts: In the inner region near the black hole the metric is approximately 
the Schwarzschild metric with a small perturbation caused by the external spacetime 
through which it is moving. In the outer region far from the black hole the metric 
is approximately the background geometry of the external spacetime with a small 
perturbation caused by the black hole. Let a length scale of the background be 71, 
and let r be some measure of distance from the black hole. Assume that yU ^ 7?. so that 
the black hole is in a context where it is meaningful to say that its mass is small. The 
inner region extends from the black hole out to r ^ 7?.. The outer region includes all 
r ^ /i. These two regions overlap in the buffer region where ^ r <^TZ. 

When we focus on the inner region in sections El and U\ the object is a black 
hole, and we find an approximation for that consists of the singular /i/r part of 
the Schwarzschild metric plus its tidal distortion caused by the background geometry. 
Equations give a straightforward approximation for /i^^. When we focus on the 

outer region we are free to think of the object as being a point particle. Matching the 
perturbed metrics in the "matching zone," within the buffer region, in sectionlHlprovides 
an approximate solution to the Einstein equations with a remainder of 0(/i^/7?.^), which 
is uniformly valid in the limit fJ^/TZ 0, everywhere outside the event horizon as is 
demonstrated in section |H1 

The motion of the object is ultimately described as being geodesic in an abstract 
metric gab + h^b^ where /if^ is the metric perturbation which would result from a point 
particle, with the singular part h^b removed. The majority of the remainder of this 
manuscript is the elucidation of the steps which lead to the calculation of the 0(;Li) 
adjustment of a small object's world line. 

3. Locally inertial coordinate systems 

A description of the metric perturbation hab near a point mass /x moving along a 
geodesic F is most convenient with coordinates in which the background geometry looks 
as flat as possible at the location of the particle. Let 7^ be a representative length 
scale of the background geometry — the smallest of the radius of curvature, the scale of 
inhomogeneities, and the time scale for changes in curvature along F. Corresponding 
to any event p, there is always a locally inertial coordinate system for which the metric 
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and the affine connection at p are those of flat Minkowski! space, rjab- The value of the 
metric and its first derivatives at p in any coordinate system are all that is required 
to determine a locally inertial system. The construction is described, for example, by 
Weinberg JHl in his equation (3.2.12). Locally inertial coordinates at p remain locally 
inertial under an inhomogeneous Lorentz transformation. In addition, if p is the origin 
of the coordinates, then any transformation of the form x^ew = x"' + \°'hcdX^x'^x'^ is also 
a locally inertial system with the origin at p. Such an O(x^) coordinate transformation 
changes the form of the metric only by O(x^) in a neighborhood of p. 

One specialization of locally inertial coordinates, which fixes the form of the 
quadratic parts of the metric at p, are Riemann normal coordinates [19^ where the 
metric takes the form 

gab = Vab - \{Racbd ' RacLbc)x''x'' + 0(xV7^='). (20) 


Any coordinate transformation of the form 

xl,^ = x'' + \\,dex'x''x''x' + 0(xV7^^) (21) 

preserves this Riemann normal form of the metric. The coordinate location of an event 
q is given in terms of a set of direction cosines, with respect to orthonormal basis vectors 
at p, and the change in affine parameter along a geodesic from p to q. Riemann normal 
coordinates are defined only in a region where the geodesies emanating from p do not 
intersect elsewhere in the region. 

Coordinates x"" = {t, x, y, z) may be found which are locally inertial along any 
geodesic F, with t measuring the proper time s on F. In these coordinates Qab = 
Vab + 0(r^/7^^), where r'^ = + y^ + = x'^Xi and the indices z, j, fc, Z, p, q run 
over the spatial coordinates x, y and z. A coordinate transformation of the form 
^new = x"" + \°'ijk{s)x^x^x'' + 0(r^/7^^) prcscrvcs these features with most components 
of the metric changing by 0(r^/7?.^). However, Qu changes only by 0{r^/TZ^) and is 
always of the simple form gu = — 1 — RtujX^x^ + 0{r^ /IZ^), where Rutj is evaluated on 
F. 



3.1. Fermi normal coordinates 

Fermi normal coordinates jED] are one specialization of locally inertial coordinates on 
a geodesic F for which the 0{r^/TZ'^) parts of the metric have a particularly appealing 
form as simple combinations of components of the Riemann tensor evaluated on F, ^] 

gab dx" dx^ = - (1 + Rtitjx'x^) dt^ - ^Rukjx'x^ dt dx'' 
+ ifki - ^Rkiijx'x^) dx'' dx' 

+ 0(rV7^^). (22) 

Li and Ni |2I] give the form of the metric in Fermi normal coordinates to higher order. 
The defining characteristics of Fermi normal coordinates are that they are orthogonal 
on F, that the spatial axes are geodesies, and that the distance from F at proper time s 
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to an event (t = is {x'^x^SijY^'^, when measured along a geodesic perpendicular to 

r. 



3.2. THZ Normal coordinates 

A second specialization of locally inertial coordinates on F, introduced by Thorne and 
Hartle |221 and extended by Zhang describe the external multipole moments, 
defined on F, of a vacuum solution of the Einstein equations. In these THZ coordinates 

Vab + Hab 

Vab + 2Hab + 3Hab + 0{ryn^), 



9ab 



(23) 



with 

2Habdx^dx' 



£,jx'x^{ &t^ + fkidx^ dx') + -ekpqB^x^x' dt dx^ 



20 

~ 21 L 

5 

+ 2l 



SijX X Xk ^ik^ 



dt dx'^ 



-r'^e ■B '^x^ 



dx* dx^' + 0(r77^^ 



(24) 



and 



+ '^ekpqB^jxPx'x^ dtdx^' + 0(rV7^^) 



(25) 



where Sijk is the fiat space Levi-Civita tensor. These coordinates are well defined up to 
the addition of arbitrary functions of 0{r^/7l'^). The external multipole moments £ij, 
Bij, Sijk, and Bijk are spatial, symmetric, tracefree (STF) tensors and are related to the 
Riemann tensor evaluated on F by 

(26) 



S - = 

B.r- 



■titj 1 

STF 



[dkR, 



and 



titj\ 



3 r „„„ „ iSTF 



(27) 
(28) 

(29) 



Bijk g dkRpqjt\ 

where ^"""^ means to take the symmetric, tracefree part with respect to the spatial indices. 
£ij and Bij are 0(1/7^^), while £ijk and Bijk are 0(1/7^^). The dot denotes differentiation 
of the multipole moment with respect to t along F. Thus £ij = 0(1/7^^) because TZ 
limits the time scale along F. All of the above external multipole moments are tracefree 
because the background geometry is assumed to be a vacuum solution of the Einstein 
equations. 

The THZ coordinates are a specialization of harmonic coordinates, and it is useful 
to define the "Gothic" form of the metric 

.,ab 



8"' = 



-99 



(30) 
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(31) 



A coordinate system is harmonic if and only if 

daH""^ = 0. (32) 

Zhang j2Sl gives an expansion of g"^ for an arbitrary solution of the vacuum Einstein 
equations in THZ coordinates, his equation (3.26). The terms of H"^^ in this expansion 
include 

-rab _ ^jjab ^ ^jjab ^ 0(rV7^^) (33) 



where 



and 



3^" 



X-' 



2 1 n r • 2 • 



5 

21 



X 



qkXpX 



(34) 



^ijk"^ X^ X 



= 0(r77^■ 

If r/7?. ^ 1 then Hab is approximately the trace reversed version of H 
H,b = Hab-lvabH'c + 0{ryn') 



ab 



(35) 



(36) 



and (j2Sl)-(|2Sl) correspond precisely to up to a remainder of O^r^/TV"). 

Zhang gives the transformation from Fermi normal coordinates to the THZ 
coordinates 



^thz — ^fn ^^*i^fn + g^jfc^fn-^fn^fn + 0(r /TZ ). 



(37) 



3.3. An application of THZ coordinates 

The scalar wave equation takes a particularly simple form in THZ coordinates, 

= (r/-^ - H''')dadhi:, (38) 

where the second equality follows from ()3ip and the last from fj32|) . After an expansion 
of the contractions on H"^^, this becomes 

v/^V^V,^ = v'^'dadk^ - md,d,^ - 2H%A)^ - H'^mij. (39) 
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An approximate solution ip with a point charge source is q/r. Direct substitution into 
(jH^ reveals just how good this approximation is. If ip is replaced by q/r on the right 
hand side, then the first term gives a (5-function, the third and fourth terms vanish 
because r is independent of t, and in the second term has no contribution because 
of the details given in ()34j) . and the 0(r^/7^^) remainder of H^^ yields a term that scales 



Note that the remainder 0(r/7?.^) is C°. From the consideration of solutions of Laplace's 
equation in fiat spacetime, it follows that a correction to q/r, of 0{r^/7l'^), would 
remove the 0{r /IZ'^) remainder on the right hand side. We conclude that q/r + 0(r^/TZ'^) 
is a solution of the scalar field wave equation for a point charge and that the error in 
the approximation of the solution by g/r is C^. In Ref. we show that q/r is the 
singular field ip^ for a scalar charge, up to a remainder of 0(r^/7^^). This was done by 
use of a Hadamard expansion of the Green's function. 

THZ coordinates provide elementary, approximate solutions to the wave equation 
with a singular source for vector and tensor fields as well [23] . 

4. Vector and tensor harmonics 

The forms of and sHab in (jUj) and might appear unfamiliar, but they 

actually consist of £ = 2 and 3 vector and tensor spherical harmonics and have a 
close relationship with those introduced by Regge and Wheeler ^3] in their analysis 
of metric perturbations of Schwarzschild black holes. This relationship is clarified with 
an example of Sij, whose Cartesian components are symmetric, tracefree, and constant. 
However, the spherical-coordinate component Srr has the angular dependence of a linear 
combination of the Y^m's for i = 2. Thus, it is convenient to define S^'^'^ = SijU'^W , where 
n* is the unit radial vector in fiat space. E^"^^ is a scalar field which carries all of the 
information contained in the constant Cartesian components of Sij and may be used to 
generate related quadrupole vector and tensor harmonics. 

For the angular components of vectors and tensors, we find it convenient to follow 
Thome's description of the pure-spin vector and tensor harmonics jSEj, which are closely 
related to the harmonic decomposition used by Regge and Wheeler ^j. For example, 
the spin-1 vector harmonics generated by the spherical harmonic function Yim are the 
even parity 



as 0(r/7^^). Thus, 




(40) 



a 



rcrJ'VbYem 



(41) 



and the odd parity 



Y„ 



a 



—reJ'VbYim, 



(42) 



where 



CTab = gab + UaUb - UaUb 



(43) 
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is the metric of a constant t,r two-sphere, and 

Cafe = etrab, with etr9^ = = sm 9 , (44) 

is the Levi-Civita tensor on the same two-sphere. Here Ua and Ua are the unit normals 
of surfaces of constant t and constant r, respectively. 

We generalize this approach: For a vector field C,a, the parts o'a^^b which are tangent 
to a two-sphere may be described by two "potentials" ^'^^ and via 

(Ta% = raa'Vkr - re^^V^r^- (45) 

The potentials and are generally functions of all of the spacetime coordinates 
and are guaranteed to exist by the invertibility of the two dimensional Laplacian on a 
two-sphere. The factors of r are included for convenience. 

The notation for a covariant vector field is condensed by defining even and odd 
parity vectors associated with the potential ^''^ 

C ^ raa'V,r (46) 

and with the potential 

^od ^ _^e„^V,r'- (47) 

The four independent components of a covariant vector in a spherically symmetric 
geometry may be written as a sum of the form 

dx"^ = 6 dt + ir dr + + if) dx^ (48) 

in terms of the four functions ^r, ^'^^ and The capital index A is used here just 
as a reminder that the vector to which it is attached is tangent to the two-sphere. The 
A index should otherwise be considered an ordinary spacetime index in the covariant 
spirit of pK|) - (p7jl . 

Similarly for a symmetric tensor field hat, the parts which are tangent to a two- 
sphere aa'^ah^hcd may be described by the trace with respect to aah and by two potentials 
/i'^^ and h"'^ via 

- r2e(/a5)'^Ve [a/Veh'"') (49) 

The potentials /i'^^ and h"'^ are generally functions of all of the spacetime coordinates 
and are guaranteed to exist by theorems involving solutions of elliptic equations on a 
two-sphere. The factors of are included for convenience. 

The notation for a covariant tensor field is condensed by defining trace-free tensors 
tangent to a two-sphere and associated with the potential /i*^^ 

hZ = rV(/a,)^V, (a/Ve/i^^) - ^rVa.a^'^Ve (a/Ve/i'=^) (50) 
and with the potential 

hit = -r'e(/a,)'^Ve (o'/Ve/^"'^) . (51) 
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The ten independent components of a symmetric covariant tensor hat in a spherically 
symmetric geometry may be written as a sum of the form 

hab dx"" dx^ = htt dt^ + 2htr dt dr + 2 {h^A + h°^) dt dx^ 
+ h„ dr' + 2 (/i^^ + hf^) dr dx^ 

+ (^-h'^'aAB + /^AB + /^ab) dx^ dx^ (52) 

in terms of the ten functions htf, ht^, hf, h°'^, h„i h'^r i ^T^i ^^"^^i ^'^^ ^'^'^ with 
the vector field, the capital indices A and B are used here just as a reminder that the 
vector or tensor to which they are attached is tangent to the two-sphere. Otherwise, 
they should be considered ordinary spacetime indices in the covariant spirit of pn|l - (jKT|l . 

The descriptions of vector and tensor potentials in (jlSjl and (jl^ on a two- 
sphere could have been written with a derivative operator involving the usual angular 
coordinates. However, this would cloud the covariant nature of the decomposition which 
is clearly revealed above. 

The description of the vector and tensor components in terms of potentials takes 
advantage of the natural symmetry of the background geometry. For example, if a 
potential is a function of r and t times a Y^m then the resulting vector or tensor field 
is the same function times the vector or tensor spherical harmonic with the same i, m 
pair. Expressions such as the perturbed Einstein tensor take a particularly simple form 
when written in terms of the potentials in place of the components. 

We assume throughout that £ is always associated with even parity vectors and 
tensors, and that B is always associated with odd parity vectors and tensors. Thus, 
and are often understood in £^ = or i3 = B°'^. A superscript in parentheses, as in 
S^'^^ = SijU^^ , denotes the multipole index i which is also the number of indices in the 
STF tensor Sij. 

With this notation, alternative forms of and (^3)) are 

2Hab dx" dx^ = - r2£(2) ^^2 ^ ^^2 ^ ^^A ^^b^ ^ 2^B^^^ dt dx^ 
7 21 ^ 

.3 , , ^3 



+ 2^B^A^ dr dx^ - ^B^^l dx^ dx^ + 0(r77^^) (53) 



and 



^3 



iHab dx" dx'' = S^^\dt^ + dr^ + aAB dx^ dx*) 

3 



3 

+ 2-S5fM^dx^ + 0(rV7^^). (54) 
9 

5. Slowly time dependent perturbations of the Schwarzschild geometry 



When a small Schwarzschild black hole of mass /i moves through a background 
spacetime, the hole's metric is perturbed by tidal forces arising from Hab in (1221), ^^d 
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the actual metric near the black hole is 



act 
Sab 



(55) 



where the quadrupole metric perturbation 2hab is a solution of the perturbed Einstein 
equations (fT^. The appropriate boundary conditions for 2hab are that it's components 
be well behaved on the future event-horizon, in a well-behaved coordinate system, and 
that 2hab iHab in the buffer region [22j, where ^ <^ r <^ IZ. The octupole metric 
perturbation s/iaf, has a similar description. 

In [Appendix G we follow Poisson's recent analysis [231231201 of a tidally distorted 
black hole, and describe the metric perturbation for r -C 7^ in ()G.6|) - ()G.11|) . An 
expansion of the metric perturbation in the buffer region for ^ <^ r <^ TZ ultimately 
provides the even parity 

2/1^1 dx'^ dx' = - ^(2) [(^ _ 2^)2 dt^ + r2 dr^ + (r^ - 2^?)aAB dx^ dx^] 

+ [2(r + /i) dt^ + 2(r + 5/i) dr^ + (2r + 5/i)(T^ij dx^ dx^l 

15r 

- 2 ^^'^'-'^;'-^^'^ + ^^'^ g(^) dtdr + 0(/^(^)/r^), (56) 
3(r — 2yLi) 

and the odd parity 



dx" dx^ 



(2(2) 



(3r + 4/i)i3 



(2) 



dtdx^ 



ijf drdx^ + 0(/i3(2)A 



(57) 



12(r - 2/i) 

which together properly match the 0{r'^/TZ'^) terms of ()24|1 or of ()53|1 : the 0(r^/7^^) 
terms are in a different gauge. In this form £^*^2) ^nd i3*^2) a^-g considered functions of t 
and f'^2) denotes the t derivative of £^'^\ Together, these provide the quadrupole metric 
perturbation up to remainders of 0{r^/TZ^) and 0{fi^ /r^TZ^). 

The approximately time independent octupole perturbation of the small black 
hole may be treated similarly. The time independent solution of E^j^^'^l^h) = which 
is well behaved on the event horizon and properly matches the 0{r^/TZ^) terms of ()25p 
or of dSl) 

2 



shabdx"" dx 



+ 1 



3 

dr^ + 



1 - 



2fi 



r 



1 - 



2 

r — 2/ir H 

5r 



(d^2 + sin^ ^ d(/.2) 



+ 2— 
9 



r 



4/x 



1-^1 dtdx 
3r 



(58) 



The part of shab proportional to Sijk or Bijk is of 0{r^/7l'^) and not required here. 

At this level of approximation, the interactions of tidal forces with a small black 
hole have no significant effect upon the motion of the hole. From the analysis of Thorne 
and Hartle [22] the dominant tidal effect upon the motion of a nonrotating object results 
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from the coupling between the external octupole moment of the geometry Sijk and the 
internal quadrupole moment of the object Ijk] the resulting force is 

fia' ~ S'.kF", (59) 

equation (1-12) of reference For a Schwarzschild black hole, Xjk must result from 
the external quadrupole moment Sjk- With dimensional analysis we conclude that this 
tidal acceleration is no larger than 

a' ~ fj.^S'jkS''" ~ /iV^^- (60) 

This acceleration is much smaller than the 0(/i/7?.^) acceleration of the self- force which 
is the focus of this manuscript. Hence, we conclude that for our purposes the tidal forces 
resulting from (|3U ]) -(|3H |) exert no significant net force on the black hole. 



6. A small black hole moving through a background geometry 
6.1. Buffer region 

In the previous section we treated the actual metric of a small black hole moving through 
an external universe as the Schwarzschild metric being perturbed by tidal forces with a 
small perturbation parameter r/TZ, 

= ^?af™ + 2hab + 3Kb + 0(rV7^^), (61) 

The metric perturbations 2^ab and ^hab are the dominant perturbations arising from the 
quadrupole and octupole tidal forces and are given in 

In the buffer region fi <^ r <^ 71 the actual metric is described equally well as the 
background metric being perturbed by the small mass ^ with a perturbation parameter 
fi/r. With THZ coordinates the background metric is 

(7°, = r]ab + 2Hab + sHab + 0{r^ /U'') (62) 

and the actual metric is 

9Tb = gib + Kb + Kl + Kl + --- (63) 

Each h'^f^ is the part of the metric perturbation which is proportional to /i". These are 
obtained by a re-expansion of the results of the previous section in terms of powers of 
the small parameter fi/r. Thus, 

Kb^oK, + 2K, + 3K, + {f^ryn'), (64) 

where 

oKbdx'' dx^ = 2^( dt^ + dr^) (65) 

is the fj,/r part of the Schwarzschild metric gab^"", 

2K, dx'^ dx' = 4/ir£(2) ^^2 _ 2^sJ) dt dx^ 

3 

+ 2^^(2) ^ 2^i3j) dr dx^ (66) 
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consists of the fir/TZ"^ and ^r'^/TZ^ parts of 2hab from and (fHTj) . and is the 
fir'^/TZ^ part of s/ia^ in (j3Hj) 



3/1;:, dx'' dx" = ^£^'^ [5 dt' + dr^ + 2r2(d^2 ^ ^^^2 ^^^2^j 

-2i^Bj)d^d.^. 
27 ^ 



(67) 



Asymptotic matching 



To add a modest amount of formahty to this analysis, we assume that the background 
metric g^f^ with a geodesic F has an expansion in terms of THZ coordinates as in 
(jnH)- We then consider a sequence of metrics gabilA which are solutions of the vacuum 
Einstein equations with a Schwarzschild black hole "centered on F" in the sense that 
near the black hole the metric is approximately described as in (jHT|) . The sequence is 
parameterized hy ji <^ TZ with gabif^) = dab- Our focus is on the behavior of gabifJ^) in 
the limit that /i — 0. This analysis falls under the purview of singular perturbation 
theory gab{.f) has an event horizon if and only if /x 7^ 0; therefore, the exact metric 
for /i = differs fundamentally from a neighboring metric obtained in the limit yU ^ 0. 

In the buffer region gab{f) is nicely illustrated in a fashion introduced by Thorne 
and Hartle [22] as a sum of elements of positive powers of the small parameters /x/r and 
r/7^, 
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■ = h^^' (68) 



^Sc^™ 2h' 3h' ih' 

where & means "and an element of the form . . ." Starting with £ = 0, the £th column 
in the tableau consists of elements which scale as {r /TZy. Starting with n = 0, the nth 
row consists of elements which scale as (/x/r)". In the ^/TZ — *• limit, every non-zero 
element in the tableau is larger than all elements below it in the same column, or to its 
right in the same row. 

The primes on the if's in the top row work around a deficiency in our notation: 
In section IT^ the prefix 2 in refers to the multipole index ^ = 2. In the tableau, 
the prefix 2 on refers to the power of the order behavior, 0(r^/7^^). While 2-f^a6 
includes not only the quadrupole parts proportional to £ij and Bij, which are 0(r^/7^^), 
but also the parts proportional to time derivatives of £ij and Bij , which are the order of 
a higher power of r/lZ. In the tableau, the time derivative terms of iHab are included 
in i+iH'^fj and columns further to the right. 

Row n is proportional to /i" and is an expansion in the external moments and in 
their time derivatives. Each element in the tableau is a finite combination of terms 
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which scale with the same power of 1/TZ, 

fi"r'-yn' ^ (^^y r' {Se & & & ■■■ & (^-2)^2) (69) 

The prefix superscript is the number of time derivatives, and ^'^'^Ei represents the even 
or odd parity i indexed STF external multipole moment differentiated with respect to 
time p times. Thus, i is the largest external multipole index that contributes to any 
element in column i or to ih'. 

At the outer edge of the buffer region, where fi/r <^ r/TZ, QabifJ') is approximately 
the background metric perturbed by fi. In this region, the top row of the tableau 
consists of the expansion of (7°^ about T in powers of r/TZ, contains no /i dependence 
and dominates the actual metric gabif^)- The sum of the top row is (7°^. 

The n = 1 row combines to give /i^^, which is the 0(/i) metric perturbation of g^/^. 
And the nth row combines to give the 0(/i") perturbation; higher order perturbation 
theory for the background geometry is necessary to determine the n > 1 rows. 

At the inner edge of the buffer region, where njr ^ f/TZ^ Qab^lj) is approximately 
the Schwarzschild geometry perturbed by background tidal forces. The £ = column 
of the tableau is simply an expansion of the Schwarzschild geometry in powers of /x/r, 
contains no TZ dependence and dominates the actual metric gabifJ^)- 

The i = 1 column, linear in r/TZ, would be a dipole perturbation of the 
Schwarzschild geometry. But there is no r/TZ term in an expansion about a geodesic. 
Consequently the top element of the i = 1 column is zero, as are all elements of this 
column. 

The top term in the i = 2 column, represents the external quadrupole tidal 
field. When this is combined with the rest of the i = 2 column the result is 2^1?,; the 
entire quadrupole perturbation of the black hole caused by tidal forces, in the time 
independent approximation. 2^1b is given explicitly as in the 0(1/7^^) terms of (fSUj) and 

Similarly, the top term in the i = 3 column, s-ff^^ represents the 0{r^/7Z^) external 
tidal field which distorts the black hole creating 3/1^,^, which is given as the 0(1/7^^) 
terms in ()5fij) - (jHH|l . Thus, the top element of each column provides a boundary condition 
for the equations which determine the resulting metric perturbation of the black hole. 
Each column also satisfies appropriate boundary conditions at the event horizon. 

The analyses for ih'^f^ up to £ = 3 are straightforward problems in linear perturbation 
theory of a Schwarzschild black hole. The nonlinearity of the Einstein equations first 
appears in the elements of the i = 4 column, which have some contributions from terms 
quadratic in the i = 2 elements. Higher order perturbation theory for a black hole is 
necessary to determine the i > 4 columns. 

The actual metric is accurately approximated by g^f™ + 2^a6 + 3^'ab r <t^7Z, and 
g^ij is an accurate approximation of gab^jj) for jj, <^ r. In the buffer region fi <^ r <^ 7Z 
these approximations are 

9at^ + 2Kb + 3K, = Vab + 2Kb + sKb + 0(/i/r) (70) 
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and 

9l = Vab + 2K, + sK, + 0(rV7^^). (71) 

A demonstration of asymptotic matching [201 requires a matching zone, within the buffer 
region, where the smallest displayed term on the right hand side, s-ff^^ = 0{r^/TZ^), is 
simultaneously much larger than both remainder terms, 0(/i/r) and 0(r^/7?.^). The 
actual metric is accurately approximated by equation (ffn|) to the "left" of the matching 
zone, by equation (f7T|) to the "right" of the matching zone, and by rjab + + s^'ab 
only within the matching zone. 

The matching zone is thus bounded by /i/r <^ r^/TZ^ on the left and by r^/TZ'^ ^ 
r^/TZ^ on the right. These may be combined into 

(/i7^=')^/^ < r < 7^, (72) 

and this fits within the buffer region because 

fi < (/i7^3)^/^ < r < 7^, /i/7^ ^ o. (73) 

This is the signature of a matched asymptotic expansion. 



7. Singular field h% 

The Einstein tensor is the sum of terms consisting of the product of various components 
of the metric and its inverse along with two derivatives. In the buffer region, where 
/i ^ r ^ 7^, an expansion of the Einstein tensor Gab[g{fJ')] in positive powers of /i/r and 
r/TZ may be represented in a tableau similar to that for gabif^) introduced in section IHl 

In the expansion of Gab[g{lJ')] the terms of every power of 1/7^ which contain no 
dependence upon /i are each zero because (/^^ is assumed to be a vacuum solution of the 
Einstein equations, Gab{g^) = 0. Similarly, all of the terms linear in /i must combine to 
yield E^bih'') = -SnTab, because /i^^ is a perturbative solution of the Einstein equations 
with Tab representing a point mass. The individual terms in gabit^) which are linear in 
jj, also form an asymptotic expansion for h'^; these are the ih'^^ terms in the n = 1 row 
of the tableau for gatifJ')- 

In sections ^ and |21 we discussed the actual metric perturbation h'^f from a point 
mass moving through an external geometry. The Hadamard form of the Green's function 
for the operator Eah{h) provides a decomposition h^^ = h^b + h^f, in a neighborhood of 
r, where Eab{hP) = —STcTab- The analysis of the Green's function yields an asymptotic 
expansion for h^f^. The remainder /i^^ is necessarily a vacuum solution of Eab{h^) = 
in a neighborhood of F where an expansion for h^b is regular. Thus, in the tableau 
for Oabit^), h^b is 0{^). However, its regular behavior in a neighborhood of F implies 
that it has no spatial dependence on a scale of 0(/i), and that it should properly be 
moved up in the tableau to be absorbed in the definition of g^b- This 0(/i) change in 
g^b would affect the /i^^ only for n > 2. Further, the actual constructions of 2h'ab and 
sh'^by resulting in equations do not appear to allow for the inclusion of any 

such regular part, except in the top row. 
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The possibility that /i^^ when promoted to the top row, might contain a dipole part 
in the i = 1 column is discussed in section |H| 



and that /i^^ is included in the top row of the tableau (jHHjl . We have verified that o^ab 
and 2h'^b (^^ Lorenz gauge) are equivalent to the first two terms in the expansion of 
h^f^ via the Hadamard form of the Green's function. Further, the ih'^l have no dipole 
i = 1 component at 0{fi) which could effect the world line F at a level of interest in a 
self-force calculation. 

In the next two sections the effect of coordinate choices on the form of /i^^ 
are discussed. First, a change in the locally-inertial coordinates appears as a gauge 
transformation of the Schwarzschild metric being perturbed by the external tidal fields. 
Second, an 0(/i) coordinate change appears as a gauge transformation of the background 
metric being perturbed by a point mass fi. 

7.1. Coordinate transformations of the locally inertial coordinates 

The convenient THZ coordinate system is used in sections El and to determine the 
leading terms oh^b, 2^^^ and 3/1^^ in an expansion of /i^j,. But, if /i^^ is to play a 
fundamental role in radiation reaction and self-force analyses then the definition of 
h^jj should certainly not be wed to any particular locally-inertial coordinate system. 

In this section we examine the change in the description of h^j^ under a change 
of locally-inertial coordinates. The next section describes how an 0(/ir^/7^^) gauge 
transformation of the perturbed Schwarzschild metric changes the form of /i^^ while 
remaining with the same locally-inertial coordinates. 

For the "inner" perturbation problem of the matched asymptotic expansions, the 
external tidal fields are considered a perturbation of the Schwarzschild geometry. From 
this perspective a change from one locally-inertial coordinate system to another appears 
as a gauge transformation of the perturbed Schwarzschild metric. 

A second locally-inertial coordinate system is defined by 



where X"'ijk is an 0(1/7^^) constant, such as in equation ()37|) which relates Fermi 
normal to THZ coordinates. For the perturbed Schwarzschild metric this appears as 
a gauge transformation with a gauge vector = X°-ijkX'^x^ + 0{r^ /TZ^). Under such a 
change in coordinates the description of /i^^ changes in two different ways: the functional 
dependence upon coordinate position changes and the components of the tensor change. 
Let the components in the y coordinate system be denoted by a prime. For a fixed 
coordinate position k'^. 




= + Wjkx'x^x^ + 0(rV7^^) 



(75) 




dx"" dx'^ 



+ 0{^ir^/n^) 



(76) 



Qya' Qyb' 
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which, when expanded out, is 

hlb' = hi, - Cdchl - 2hl^A)C + 0(/irV7^3). (77) 

The left hand side is evaluated at y'^ = /t^ and the right hand side at x'^ = k,'^. In terms 
of the Lie derivative £, the descriptions of the single tensor field /i^^ in two different 
locally-inertial coordinate systems are related by 

/if,,, = /if, - £^ hi + 0(/irV7^3). (78) 

Now, h^f^ = 0/1^6 + 2/iafe + 0(/ir^/7^^), as in and o/^af, = ^{l^/'^) ^'^Y locally-inertial 
coordinates. Thus, the change in hl^^ is most naturally assigned to 2/iafe) 

2hT^ = 2hT - oKb + 0(/.rV7^3). (79) 
This description of the change in the 2h'^b part of /i^j, is consistent with the related 
gauge transformation of the i = 2 metric perturbation, 2/^06 = 2-/^afe + 2/^afe + 0(/i^/7^^), 
of the Schwarzschild geometry 

2Kr = 2hf,^ - £^9lt^ + 0(rV7^3). (80) 
The leading terms of this for large r are 

2HT + = 2Ht^ + - £^{vat + oKt) + 0{ryn\ /irV7^^). (81) 

These are naturally apportioned as 

,H:r = 2h:^ - £^vat + o{ryn') 
= 2^ 

and 

2Kr = 2Kf - oKb + 0(/.rV7^^). (83) 

A comparison of (f7^ and (jH^ reveals the consistency of the description of /i^, as 
a single tensor field, which in any normal coordinate system is approximated by 
oKb + 2Kb + A + 0(/xrV7^^) for /x « r « 7^. 

An 0[r^/TZ^) transformation changes 3/1^, in a similar way. 

7.2. Transformation of h^^^ to the Lorenz gauge 

The convenient Regge- Wheeler gauge was used, with the THZ coordinates, in sections 
Elandinito determine the leading terms oh'^^^, 2/1^^ and 3/1^^ in an expansion of /i^,. But, 
if /i^, is to play a fundamental role in radiation reaction and self-force analyses then the 
definition of /i^, should certainly not be wed to any particular gauge choice. 

This section gives an example of an 0{fir^ /TZ^) gauge transformation, for i = and 
2, of the perturbed background metric g^ which changes the form of h^ while remaining 
with the same locally-inertial coordinates. The previous subsection describes how the 
description of /i^j, changes under a change of locally-inertial coordinates. 

h^ is given above in in the Regge- Wheeler gauge. To transform the 2h'^^ 

part of this into the Lorenz gauge, the gauge vector is 



,H:^ - 2V(„6) + 0(rV7^3) (82) 
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The Lorenz gauge has 

2K» = 2KU^W) - - V,ea, (85) 

where the metric being perturbed is (7°^, and Va is its covariant derivative operator. 
This resuhs in 

2Kli\z) dx" dx' = ^ [(1 + r''£^'^) dt' + (1 - 3r2^(2)) dr' + (1 - r^£^^^)aAB dx^ dx^] 
- 4^r^^^^ dr dx^ - dx^ dx^ 

+ 2|^i3J^d^dx^ + 0(/irV7^3). (86) 

For completeness, the trace of 2^ab(lz) is 

(^af _ 2/7-^)2/i^;(lz) = 4/x/r + 0(/xrV7^3), (87) 
and the trace-reversed 2h'^b = "i^ab ~ l9lb9o^2K'd is 

2/i^;(lz) dx^dx' = ^(1 + r2^(2)) ^^2 _ 4^^^(2) ^^2 

- 4/ir4^) dr dx^ - 2firS^^l dx^ dx^ 

-2^Sfdtdx^, (88) 
3 

which satisfies the Lorenz gauge condition. 

V^^::;(lz) = 0{^^r/n'). (89) 

Equation (jHUj) gives 2^0^ in the Lorenz gauge with THZ coordinates. From 
the perspective of the background metric g^f^, a change from THZ to Fermi normal 
coordinates, as described in the previous section, would preserve the covariant condition 
(jHn|l for the Lorenz gauge and provide /^^^(Iz) in Fermi normal coordinates. 

8. Regular field /i^ 

In a self-force application, it is first required to find the actual metric perturbation h^f' 
for a point mass fi moving along a geodesic F of the background spacetime g^f^. In many 
cases h'^f will be the retarded metric perturbation. However, we prefer to leave the 
choice of boundary conditions general. 

From the expansion of (7°^ about F, as in El the first few terms of an asymptotic 
expansion for h^f^ is determined as in El The regular remainder is defined by 

hi ^ - hi (90) 

in a neighborhood of F where Eab{h^) = 0. h^^, does not change over an 0{n) length 
scale, so it is natural to combine /i^j, with g'^^ in the top row of the tableau of IHl Then 
the condition that the dipole term of the top row is zero is equivalent to the condition 
that F is actually a geodesic of g^f^ + h^f^. 

From a different perspective, if /i^^ is left in the /x^ row, and if F is not a geodesic of 
9ab + ^ab; then necessarily has a dipole part in its expansion about F. This implies 
that the gravitational field of /i is not centered upon F. The act of adjusting F to 
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remove the dipole field and to accurately track the center of the gravitational field of 
fi is equivalent to requiring that F be a geodesic of g^f^ + h^^. This act of adjustment 
is also equivalent to performing a gauge transformation to the perturbed Schwarzschild 
geometry that removes the dipole field. 

Thus the consistency of the matched asymptotic expansions implies that, indeed, 
the 0(/x) correction to geodesic motion for an infinitesimal black hole has the motion 
being geodesic in g^^ + h^^. 

In an actual calculation, an exact expression for h^^^ is usually not available. It 
is only necessary that (7°^ + h^^^ be in order to determine a geodesic, and this 
requirement can be met as long as the approximation for /i^^ includes at least 
q/i^' + 2h^' ■ The next term is 3/1^' = 0(yur^/7^^), and its derivative necessarily vanishes 
on r where r = 0. Thus, calculations of the self-force will be successful as long as the 
monopole and quadrupole terms of the asymptotic expansion for /i^^ are included in the 
evaluation of h^^^ via (jHOJ. Nevertheless, including the higher order terms of /i^^, results 
in the approximation for /i^^ being more differentiable. In a calculation, usually h'^^ 
is determined as a sum over modes with /i^^ being decomposed in terms of the same 
modes. In determining the self-force, the more differentiable the /i^^ is, the more rapidly 
the sum over modes converges. The use of higher order terms in an approximation for 
h^^fj can have dramatic effects on the ultimate accuracy of a self-force calculation 

If we have the actual 0(/i) metric perturbation for a point mass, then the 
asymptotic matching provides an approximation for the geometry of a small black hole 
moving in the external geometry gabifJ') in the limit that n/TZ 0. The approximation 
extends throughout the entire external spacetime down to the event horizon. Further, 
the approximation is revealed to be uniformly valid by the concise description of the 
matched geometry as 

9ab{fi) = {9I + Kb) + {glt^ + 2K, + ,h'j 

- {r^ab + 2HU + sKb + oKl + 2KI + shO + 0{f,'/n'), /i/7^ 0. (91) 

The combination g^f^ + h^f includes all terms in the top two rows of the tableau 
but extends outside the buffer region to include the entire external spacetime. The 
combination g^^^"" + 2^af) + s^afe the left four columns in the tableau. The remaining 
terms keep the entire expression from double- counting the elements in the upper left 
corner. The dominant term from the tableau which is not included here is 0(/i^r^/7^^) 
which gives the 0{fi^/TZ^) remainder for this uniformly valid approximation for the 
matched metric in the limit fJ^/TZ 0. 

9. Gauge issues 

9.1. Gauge transformations 

In perturbation analyses of general relativity jnH E21 ISH] , one considers the difference 
in the actual metric g^f' of an interesting, perturbed spacetime and the abstract metric 
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g^f^ of some given, background spacetime. The difference 

Kh = 9Th - 9l (92) 

is assumed to be infinitesimal, say 0{h). Typically, one determines a set of linear 
equations which govern hah by expanding the Einstein equations through 0{h). The 
results are often used to resolve interesting issues concerning the stability of the 
background, or the propagation and emission of gravitational waves by a perturbing 
source. 

However, is ambiguous: The metrics g^l^ and g^^ are given on different 

manifolds. For a given event on one manifold at which corresponding event on the other 
manifold is the subtraction to be performed? Usually a coordinate system common to 
both spacetimes induces an implicit mapping between the manifolds and defines the 
subtraction. Yet, the presence of the perturbation allows ambiguity. An infinitesimal 
coordinate transformation of the perturbed spacetime 

x'" = x^ + where C = 0{h), (93) 

not only changes the components of a tensor at 0(/i), in the usual way, but also changes 
the mapping between the two manifolds in (jHS))- After the transformation (jHSI), 

f^ab - [9cd + f^cd ) Q^,a [9ab + 4 g^^J- [M) 

The ^'^ in the last term accounts for the 0{h) change in the event of the background 
used in the subtraction. After an expansion, this provides a new description of hab 

Lnew /,old ^0 n de ,Mab 

^ab - Kb 9cbg^ 9cbg^ 4 ^ 

= hf^ - £^gl, = - 2V^a^h) (95) 

through 0{h); the symbol £ represents the Lie derivative and is the covariant 
derivative compatible with g^^^. 

The action of such an infinitesimal coordinate transformation is called a gauge 
transformation and does not change the actual perturbed manifold, but it does change 
the coordinate description of the perturbed manifold. 

A similar circumstance holds with general coordinate transformations. A change in 
coordinate system creates a change in description. But, general covariance dictates that 
actual physical measurements must be describable in a manner which is invariant under 
a change in coordinates. Thus, one usually describes physically interesting quantities 
strictly in terms of geometrical scalars which, by nature, are coordinate independent. 

In a perturbation analysis any physically interesting result ought to be describable 
in a manner which is gauge invariant. 

9.2. Gauge invariant quantities 

Gauge invariant quantities appear to fall into a few different categories. 
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The change in any geometrical quantity under a gauge transformation is determined 
by the Lie derivative of that same quantity on the background manifold. This is 
demonstrated for the gauge transformation of a metric perturbation in We also 

used this fact to describe the change in /i^^ under gauge transformations in sections 
17.11 and 17.21 Thus, if a geometrical quantity vanishes in the background, but not in 
the perturbed metric, then it will be gauge invariant. Examples include the Newman- 
Penrose scalars \E'o and which vanish for the Kerr metric. In the perturbed Kerr 
metric \E'o and are non zero, gauge invariant and the basis for perturbation analyses 
of rotating black holes. A second example has the background metric being a vacuum 
solution of the Einstein equations, so its Ricci tensor Rab vanishes. The Ricci tensor 
of a perturbation of this metric is then unchanged by a gauge transformation. This is 
directly demonstrated in [Appendix D| 

Some quantities which are associated with a symmetry of the perturbed geometry 
are gauge invariant. For example a geodesic of a perturbed Schwarzschild metric, 
where the perturbation is axisymmetric with Killing field k"', has a constant of motion 
k'^u'^idab + ^ab) which is gauge invariant. 

Another symmetry example involves the Schwarzschild geometry with an arbitrary 
perturbation. It is a fact that a gauge transformation can always be found, such that 
the resulting hab has the components hgg, h^^ and /i^^ all equal to zero. In this gauge, 
the surfaces of constant r and t are geometrical two-spheres, even while the manifold as 
whole has no symmetry. The area of each two-sphere can be used to define a radial scalar 
field R which is constant on each of these two-spheres. This scalar field on the perturbed 
Schwarzschild manifold is independent of gauge. However, its coordinate description in 
terms of the usual t, r, 6, cj) coordinates does change under a gauge transformation. We 
find a use for this gauge invariant scalar field in section 110.31 

Quantities which are carefully described by a physical measurement are gauge 
invariant. For example, the acceleration of a world line could be measured with masses 
and springs by an observer moving along a world line in a perturbed geometry. The 
magnitude of the acceleration is a scalar and is gauge invariant. If the world line has 
zero acceleration, then it is a geodesic. Therefore, a geodesic of a perturbed metric 
remains a geodesic under a gauge transformation even while its coordinate description 
changes by 0(/i). 

The mass and angular momentum are other gauge invariant quantities which might 
be measured by distant observers in an asymptotically flat spacetime. A small mass 
orbiting a larger black hole perturbs the black hole metric and emits gravitational waves. 
The gravitational waveform measured at a large distance is also gauge invariant. 

9.3. Gauge transformations and the gravitational self-force 

We understand that a point mass moves along a world line of a background metric g^^^ 
and causes a metric perturbation h^, which may be decomposed into h^^ and /i^^. The 
gravitational self-force makes the world line be a geodesic of (7°^, + /i^j,. This world line 
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is equivalently described in terms of the background metric and its perturbation by 

u'V.u^ = -{gf + u^u'Xu\V.hl - ^V,hf,) = 0{f,/n') (96) 

where the covariant derivative and normahzation of are compatible with (7°^. 

Given this world line, let be a differentiable vector field which is equal, on the 
world line, to the 0(/i) displacement back to the geodesic of (7°^ along which the particle 
would move in the absence of /z^^; otherwise is arbitrary. Such a generates a gauge 
transformation for which the right hand side of ()9fij) is zero when evaluated with the 
new hab With the new hab there is no "gravitational self- force", and the coordinate 
description of the world line is identical to the coordinate description of a geodesic of 
g^f^. With or without the gauge transformation, an observer moving along this world 
line would measure no acceleration and would conclude that the world line is a geodesic 
of the perturbed metric. 

This example shows that simple knowledge of the gravitational self-force, as defined 
in terms of the right hand side of (jHEj), is not a complete description of any physically 
interesting quantity. 

With this same example, after a time T the gauge vector ^ ~ T^m ~ T'^fi/TZ'^, and 
as long as T < 7^ the gauge vector ^'^ remains small. However, when T ~ lZ^/lZJJ^, 
which is much larger than the dynamical timescale TZ, the gauge vector ^ ~ 7?. and 
can no longer be considered small. Thus, a gauge choice which cancels the coordinate 
description of the self-force necessarily fails after a sufficiently long time. Mino [31] takes 
advantage of this fact in his proposal to find the cumulative, gravitational self-force effect 
on the Carter constant for eccentric orbits around a rotating black hole. 

10. An example: self force on circular orbits of the Schwarzschild metric 

The introduction described the Newtonian problem of a small mass /i in a circular orbit 
of radius R about a much larger mass M. The analysis results in the usual 0(/i/M) 
reduced mass effect on the orbital frequency f2 given in Q. Reference j2] has a thorough 
introduction to the mechanics of this self-force calculation and gives a detailed discussion 
of this elementary problem using the same language and style which is common for the 
relativistic gravitational self-force. This includes elementary expressions for the S and R- 
fields of the Newtonian gravitational potential with descriptions of their decompositions 
in terms of spherical harmonics. 

The extension of this Newtonian problem to general relativity is perhaps the 
simplest, interesting example of the relativistic gravitational self- force. Thus, we focus 
on a small mass /i in a circular geodesic about a Schwarzschild black hole of mass M, and 
we describe each of the steps necessary to obtain physically interesting results related 
to the gravitational self-force. 
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10.1. Mode sum analysis 

Metric perturbations of Schwarzschild have been thoroughly studied since Regge and 
Wheeler [TnilII|. Both Tab and hab are fourier analyzed in time, with frequency u, and 
decomposed in terms of tensor spherical harmonics, with multipole indices i and m. 
Linear combinations of the components of /i^™'"^ satisfy elementary ordinary differential 
equations which are easily numerically integrated. With the periodicity of a circular 
orbit, only a discrete set frequencies Um = —mfl appear. 

We assume, then, that h^JI^''^{r) can be determined for any i and m. The sum of 
these over all i and m then constitutes h^f, and this sum will be divergent if evaluated 
at the location of /x. 

The next task is to determine h^f^. The THZ coordinates, including 0{r^/TZ^) terms, 
for a circular orbit of Schwarzschild are given in reference jSlj. Equations (pTK|) - ()(i7|l give 
an approximation for /i^j, in THZ coordinates with a remainder of 0(/ir^/7^^). 

We follow the mode-sum regularization procedure pioneered by Barack and Ori 
[Sni IH^ inn IMI and Mino, Sasaki and Tanaka [HHl HO] and followed up by others 
[211 im 112] ■ In this procedure, the multipole moments of the S-field are calculated 
and referred to as regularization parameters. The sum of these moments diverges when 
evaluated at the location of /i, but each individual moment is finite. Importantly, the 
S-field has been constructed to have precisely the same singularity structure at the 
particle as the actual field has. Thus the difference in these moments gives a multipole 
decomposition of the regular R-field. Schematically, this procedure gives 



for the regular field. 

We note that the sum over modes of a decomposition of a C°° function converges 
faster than any power of i. And, the less the differentiability of the function then the 
slower the convergence of its mode sum. Exact values for would then give rapid 

convergence of the sum yielding a C°° representation of /i^^. However, the approximation 
for /i^j in ()65 p - (jU7j) has an 0(/ir^/7?.^) remainder which is necessarily only C^. This 
immediately puts a limitation on the rate of convergence of any mode sum for /i^^. 
Further, /i^^ is only defined in a neighborhood of yU. Whereas a decomposition in terms 
of spherical harmonics requires a field defined over an entire two-sphere. It is important 
that the extension of /i^^ over the two-sphere is C°° everywhere, except at /i, to insure 
rapid convergence of the mode sum. This ambiguity for /i^j,, away from /i, highlights 
an important fact: the value of any individual multiple moment fi^^^^'^'^ or f^^^^^'^^ ig 
inherently ill defined. Only a sum over modes, such as in ()97|) . might have physical 
meaning. 

10.2. Gauge issues 

A vexing difficulty with equation (jU7|) revolves around gauge transformations. What 
assurance do we have that the singularity structure of h^^^ truly matches the singularity 




(97) 
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structure of /i^^*? For example /i^^ is often described in the Lorenz gauge which is well 
behaved by most standards, whereas h^f' is most easily calculated in the Regge- Wheeler 
gauge which often entails discontinuities in components of hab- 
A gauge transformation does not change the relationship 

Eabih^"') = Eatih'') + Eatih"") = -SvrT^b. (98) 

But it also does not dictate how to apportion a gauge transformation for h'^f between 
h^ij and /i^^. In a neighborhood of the particle /i^^ is known to be a solution of 
Eab{h^) = 0, but a gauge transformation generates a homogeneous solution — 2V(a^b) to 
the same equation, thus /i^^ can determined only up to a gauge transformation. Even a 
distribution- valued gauge transformation might be allowed allowed because [Appendix F 
shows that Eab(VO = 0, in a distributional sense, even in that extreme case. Thus it 
is expected that h^j^ calculated from ^7\i might have a non-differentiable part resulting 
from a singular gauge difference between and /i^^. 

My personal perspective on this situation is reassuring, at least to me, but certainly 
not rigorous. I have considered about a half-dozen different gravitational self-force 
problems involving a small point mass orbiting a much larger black hole. In each problem 
the goal is the calculation of an interesting, well-defined gauge invariant quantity. For 
each of these, the natural formulation of the problem shows that there are ways to 
define and to calculate the relevant quantities which are not deterred by a difference in 
gauge between /i^^ and h^f', even if the difference involves a distribution- valued gauge 
vector. It appears as though a particularly odious gauge choice might exist for a specific 
problem, which might interfere with a calculation. However, none of a wide variety of 
natural choices for a gauge have this difficulty for the problems that I have examined. 
Specifically, the example in this section appears to avoid any difficult gauge issue. 



10.3. Geodesies of the perturbed Schwarzschild metric 

A particle of mass /i in a circular orbit about a black hole perturbs the Schwarzschild 
metric by h'^^ = 0(/i). The circumstances dictate boundary conditions with no 
gravitational radiation incoming from infinity or outgoing from the event horizon. 

The dynamical timescale for a close orbit is 0(M) and much shorter than the 
dynamical timescale 0(M^//x) for radiation reaction to have a significant effect upon 
the orbit. Thus the particle will orbit many times before its orbital frequency changes 
appreciably. Under these conditions, the perturbed metric appears unchanging in a 
coordinate system that rotates with the particle. Thus, for times much less than the 
radiation reaction time, there is a Killing vector A;", 

£k{9'ab + h:^^) = 0, (99) 

whose components in the usual Schwarzschild coordinates are 
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Let an observer on the particle be equipped with a flashhght which he holds pointing 
in the plane of the orbit at a fixed orientation with respect to the tangent to the orbit. 
In other words the orientation of the flashlight is Lie derived by k"" and the beam of 
light sweeps around the equatorial plane once every orbit. A distant observer, in the 
equatorial plane, measures the time AT between the arrival of two flashes of light from 
the particle and concludes that Q = 27r/AT for the particle when the light was emitted. 
This operational measure of fl is independent of any gauge choice made for h'^f. 

The components of the Killing vector fc" in (jlOOj) are actually only correct in a 
particular gauge for which both £kgab — ■^kh^b — 0; individually. Under a gauge 
transformation the coordinate description of k"" changes by 0(yu) with A A;" = — A 
choice for for which f^k"' is not zero is allowed but it would be very inconvenient and 
would result in a gauge for which both £fc5'°b = 0{h) and £kh^b — 0(^)5 even though 
dnni) would still hold. In principle, the geodesies of the light rays could be computed 
from the particle out to the distant observer in this inconvenient gauge and the orbital 
period could still be determined. But in practice, this task would be horrendous. 

In the convenient gauge, with £kh^b ~ ^1 calculation is much easier. From 
the symmetry, it is clear that the change in Schwarzschild coordinate time between the 
reception of two light flashes at the observer is the same as the change in Schwarzschild 
coordinate time at the emission of these flashes. Thus, AT measured operationally by 
a distant observer is equal to the AT at the particle for one complete orbit, as long as 
a gauge is used which respects the inherent symmetry of the example. 

We next derive an expression for Vt in which is explicitly gauge invariant 

for any transformation which respects the symmetry of the example. This includes the 
possibility of a singular gauge transformation of the type that would transform from 
the Regge- Wheeler gauge to the Lorenz gauge. 

We let the particle /i move along a geodesic of the perturbed Schwarzschild 
geometry, gab + hab, where hab is the regular remainder h^^^ for /i in a circular orbit 
in the equatorial plane. The geodesic equation for the four-velocity of /i is 

^ = ^'"'a^^^^^ + ^^^) 
The perturbation breaks the symmetries of the Schwarzschild geometry, and there 
is no naturally deflned energy or angular momentum for the particle. However we let 
R{s) be the value of r for the particle, and we deflne speciflc components of Ua by 

ut = -E, u^ = J = R, and = {) (102) 

where ' denotes a derivative with respect to s. E and J are reminiscent of the particle's 
energy and angular momentum per unit rest mass. 

The components of the geodesic equation ()1U1|) are 

^ = -i.V^ (103) 
ds 2 dt ^ ' 

ds 2 d(h ^ ^ 
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We are interested in the case when the orbit is nearly circular with R only resulting 
from the effects of energy and angular momentum loss. In this case, E and J are 0{h), 
and we look for the additional condition that R is also 0{h) to describe the slow inspiral 
of /i. All of the following equations in this section are assumed to be correct through 
0(/i), unless otherwise noted. 

The normalization of is a first integral of the geodesic equation, and with the 
assumption that R = 0{h) this is 

- W'u'igab + hab) = 1 = ^_2M/j^ - ^ + W'u'hab. (106) 

While neither d/dt nor d/d(j) is a Killing vector of Qab + hab, the combination, 
k°- = d/dt + Qd/d(f) is a Killing vector in a preferred gauge, and is tangent to a 
trajectory of this Killing vector, up to 0{h). Thus, at a circular orbit u°'dahbc = 0{h?) 
in Schwarzschild coordinates. 

A description of the quasi- circular orbits is obtained from fll05|) and ()106|) by setting 
R to zero. The results are 

= - - iRu^u'drh^b) (107) 

and 

, MR^ , n b, . R^(R-2M) . ^ , , , 

' = R^'-' - ^" " - 2(fi-3M) '"' '-"°'- 
Also the angular velocity, Q, of a circular orbit as measured at infinity is 

= (d0/dt)2 = {u*/y'f = M/R' - ^~^^ u-u%hab. (109) 

Finally, 

{E - njf = (1 - 3M/i?)(l - u^v^hab + Ru^u^drhab/^). (110) 

These equations give E, J, VL and E — QJ for a circular orbit in terms of the radius 
of the orbit R and the metric perturbation hab- We can consider the effect on these 
expressions of a gauge transformation which preserves the d/dt + Vld/dcf) symmetry 
of the problem. The analysis uses descriptions of gauge transformations found, for 
example, in references ^H] and Here we present only the results. 

The orbital frequency Vl and E — QJ are both invariant under a gauge 
transformation, while E and J are not. However, both dE/ds and dJ/ds in ()103p 
and p()4|l are gauge invariant. This latter result might have been anticipated by using 
an operational definition of energy and angular momentum loss as measured by a distant 
observer, and by finding a relationship which joins the right hand sides of ()103|) and 
fll04j) with the matching conditions at the particle for the differential equations which 
describe the metric perturbation. This relationship is straightforward but quite tedious 
to demonstrate directly. 
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The energy and the angular momentum measured by a distant observer are gauge 
invariant. At zeroth order in the perturbation the energy is just M, the mass of the 
black hole. The i = metric perturbation gives the 0(/i) contribution to the energy 
and this is just fiE, which is an 0(/i) quantity independent of gauge but not relying 
upon the 0(/i) terms in ()l()7j) . The contribution of the gravitational self-force to the 
energy measured at a large distance shows up only at second order in /x; to calculate 
this effect requires going to second order perturbation analysis. Similar statements hold 
for the angular momentum measured at a large distance and J. 

For a circular orbit the radius, R, and both E and J all depend upon the choice 
of gauge for hab- However, Q is defined in terms of a measurement made at infinity, 
and E — QJ is the contraction of with the Killing vector hence, these latter 
two quantities are independent of the gauge, and this has been demonstrated explicitly 
allowing for distribution valued gauge transformations. 

The gauge invariance of Q has an interesting twist. While Q is gauge invariant, 
the Schwarzschild radius of the orbit is not. A typical gauge vector has a radial 
component which changes the coordinate description R of the orbit. This affects Q 
through the M/R^ term in p()9j) . This radial component of C,"" also changes the u°'u^drhab 
in a manner that leaves the right hand side of p09p unchanged. Equation p09p gives 
the same result whether evaluated in a limit from outside the orbit or inside, in the 
event that hab is not differentiable at the orbit; this result follows from analysis of the 
jump conditions on hab at the orbit. 

11. Future prospects 

Within the next year or two important applications of gravitational self-force analyses 
will be viable. 

For some time, it has been possible to calculate energy and angular momentum 
loss by a small mass in an equatorial orbit about a Kerr black hole using the Teukolsky 

formalism, which involves the Newman- Penrose 0^ scalars ipQ and ip4. For orbits 
off the equatorial plane this is not good enough. For gravitational waveform prediction, 
it is also necessary to calculate the dissipative change in the third "constant" of the 
motion, the Carter constant C, due to gravitational radiation. Energy and angular 
momentum loss may be determined by finding the flux at a large distance in a gauge 
invariant manner. There is no corresponding "Carter constant flux." However, Lousto 
and Whiting jlSl HE] describe progress in determining metric perturbations from ipo or 
i/j^. And Mino jSl] has proposed a method for determining dC/dt which depends upon 
these metric perturbations. 

Self-force calculations in the Schwarzschild geometry are much easier, and progress 
is likely to be rapid both in connecting results with post-Newtonian analyses and in 
tracking the phase of gravitational radiation in an extreme mass-ratio binary. 
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11.1. Evolution of the phase during quasi- circular inspiral 

One application of self-force analysis is to track of the phase of the gravitational wave 
from a small object while it orbits a large black hole many times. For this task, Burko 
[TTj has emphasized the necessity of using higher order perturbation theory to calculate 
properly the effects of the conservative part of the self-force. Here, we follow his lead, 
and estimate the number of orbits which can be tracked by use of analysis with different 
levels of sophistication. 

For definiteness, we assume that a small mass ^ is undergoing slow, quasi-stationary 
inspiral about a Schwarzschild black hole of mass M and that the orbit is relativistic so 
that M gives the dynamical time scale. A gauge-invariant E of the orbiting particle is 
defined in terms of the mass as measured at infinity, 

fiE = Mo^- M. (Ill) 

If we know Q{E) and also dE/dt, then the assumption of quasi- circular inspiral provides 

0. The phase of the orbit is then 



dt dE dt 
Let Qo be the orbital frequency at t 

(j){t) = [ n{E{t)) dt 



(113) 



n^ + t 



dQ dE 
dE~dt 



+ 



dt 



after a Taylor expansion. 

At the lowest level of approximation Q and E are given by the geodesic equation in 
the Schwarzschild metric. The solution of the first-order metric monopole perturbation 
problem, via Regge- Wheeler jT^I analysis, gives 

E ^ i^ist = — nt(circular orbit) (114) 

First order analysis, of the sort that was available in the 1970 's, also allows for the 
determination of {dE/dt)ist. To obtain new information |17] regarding E, Q and dE/dt, 
requires second-order perturbation analysis, which presupposes the solution of the first- 
order self-force problem. Second and higher order analysis would provide 



dQ 
d^ 



dn dE 
dE^ 



[I + A2nd + 



;ii5) 



1st 



where A2nd = 0{fi/7V). With second or higher order analysis, the phase is 



tn^ + -t^ 



t 



dQ dE 
dE dt 
dn dE' 
lEltt 



[1 + Asnd + 



J 1st 



dt 



[1 + Aond + ■ 



(116) 



1st 



after integration. 

Consider the size of the contribution to the phase of the different terms of 



-t' 



dn dE 
'dEltt 



[1 + A2nd + 



=:t' 



1st 



1 + 



M 



2nd 



(117) 
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If radiation reaction is not included in the analysis, then none of this term, of order 
^t^/i/M^, is included. This would lead to a phase error of one full cycle after a time of 
order tdp = M^JM/ ^, which is known as the de-phasing timescale. 

If only the first-order radiation reaction term is included, then the \t^iJ? jM^ term is 
not included and leads to a phase error of one full cycle after a time of order t„ = M"^ / jj,. 
This is the radiation reaction timescale. 

If second-order radiation reaction is also included, then the ■ ■ ■ terms of order 
^t^jj?" /M^ are not included and create a phase error of one full cycle after a time of 
order t2nd = (M^/yu)A/M/yu). This is the second-order timescale. 

These same results are restated by noting that geodesic motion loses one cycle 
of phase information after order ^JM/ /i orbits. First order perturbation theory loses 
one cycle of phase information after order M/ fi orbits. And second order perturbation 
theory loses one cycle of phase information after order (M^//i)^M//i orbits. 

These estimates describe the difficulty involved in tracking the phase of an orbit 
over an increasing number of orbits. 

11.2. Connection with post- Newtonian analyses 

An effort is now underway to find the effects of the gravitational self-force on a number 
of parameters related to orbits in the Schwarzschild metric. The first interesting results 
will be the orbital frequency Q and the rate of precession of the perihelion for a slightly 
eccentric orbit. Other parameters which can be calculated with self- force analysis for 
circular orbits are E — QJ and a gauge invariant measure of the radius of the orbit (see 
section 19. 2j) . 

First order perturbation theory coupled with self-force analysis will provide the 
0(/i/M) effect on the innermost stable circular orbit (ISCO), as well as the effect on 
the angular frequency of the ISCO. Currently, there is no firm prediction as to whether 
the self-force moves the ISCO in or out. Some recent scalar-field self-force results |1H] 
show that the ISCO moves in and the frequency of the ISCO increases; but there is no 
clear generalization of this result to gravitation. 

More interesting quantities will be available with second order perturbation 
calculations, which now appear feasible. These include the energy, angular momentum, 
and the rate of radiative loss of these quantities. Eventually, second order gravitational 
wave-forms will be calculated. 

One early goal of self-force analysis is to make contact with post-Newtonian results. 
To do so requires that the quantities being calculated via perturbation analysis match 
up precisely with those from post-Newtonian analysis. 

Post-Newtonian analyses are most reliable with slow speeds and weak gravitational 
fields, and they easily accommodate comparable masses in a binary. Perturbation 
analyses are most reliable with an extreme mass ratio, but they accommodate fast 
speeds and strong fields. For, say, a 3Mq black hole orbiting a 20Mq black hole near its 
innermost orbit, the mass ratio is not very extreme, the speeds are not very slow and 
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the fields are not very weak. Nevertheless, for this situation both post-Newtonian and 
perturbation methods will be able to estimate properties of the system. A comparison 
of these estimates will certainly highlight the strong and weak aspects of each approach. 
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Appendix A. Perturbed Bianchi identity 

The Bianchi identity is 

'^cRdea + eRcda + ^ dReca = 0. (A.l) 

Contraction on c and b implies that 

"^bRdea' = (A.2) 

for a vacuum solution of the Einstein equations. This result is used often in the 
derivations of identities involving Eab{h). 

The definition of the operator Eab for a vacuum spacetime is 

2Eab{h) = V^hab + VaVfe/l - 2V(aV^/l,)e + ^Ra^^Kd + QabiW^h^ - V'/l), (A.3) 

SO that 

+ 2{y''Ra\'')Kd + 2R'\''yahcd + Vby^^^h^d - W.W'Wch. 

+ R°''^b'^^chad + 2R°''^b'^Vahcd 

= 0. (A.4) 

The second equality follows after use of the Ricci identity to interchange the order of 
derivatives on the first, second and last terms as well as repeated uses of Rab = and 
()A.2|1 for vacuum spacetimes. The final result follows after use of the Ricci identity 
on the first two terms and on the second two terms of the second equality, and the 
application of symmetries of the Riemann tensor on the remainder. 

If hab is not then V^Eabih) = in a distributional sense. To show this, choose 
an arbitrary, smooth test vector field A" with compact support. Consider the integral 
of X'^'V°'Eab{h) over a sufficiently large region. Integrate by parts once and discard the 
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surface term. Next use ()E.4|) and discard the surface terms to obtain an integral of 
h^'^EabiyX). This integral is zero from ()D.1|) . These steps also provide an alternative 
derivation of ()A.4jl in the event that hab is, in fact, C^. 

Appendix B. Formal nth order perturbation analysis 

In general perturbation analysis, let the gab of (fT^ be an exact solution to the vacuum 
Einstein equations, (7°^, and iteratively define 

iJab — iJab ^ "-ab \^-^) 

where 

/^i? = 0(/i"). (B.2) 
Assume that we are given g'^~^^ and T^^^^ = 0(/i), with 

Glr'^-8vrTi:) = 0(/i"). (B.3) 
If h^^^ is a solution of ()B.3|) from 

Eabih^'""^) = Glr^) - SttTI") + 0(/i"+^), (B.4) 
then it follows from the definition of the operator Eab{h) in (jl4j) that 

G'i?-8vrTi:) = 0(/i'^+^), (B.5) 

and h^^fj is an 0(/i") improvement to the approximate solution to the Einstein equations. 
The Bianchi identity implies that 

y^'Eabi^h) = (B.6) 

for any symmetric tensor field hab, as shown in [Appendix A It is also shown that if 



hab is not then ()B.6|) holds in a distributional sense. Thus an integrability condition 
of (lR4l) is that 

V"(Glr')-87rTi")) = 0(^-+i). (B.7) 
Note, however, that 

+ KciGf^' - SttT^")^) - K.iG^r'^'^ - SvtTW'^), (B.8) 

where 'V'^^_-^^ is the derivative operator of g^l ^\ and F^^ is the connection relating the 
derivative operators V" and V^„_^j. The Bianchi identity implies that 

V(„-i)Glr'^ = 0, (B.9) 
and the terms in ()B.8|1 involving F^^ are order /i""*"^ because of ()B.3|1 and the fact that 
F^^ = 0(/x). Thus, the approximate vanishing of the right hand side of ()B.8|1 is the 
integrability condition for ()B.4|) . 

V(„_i)ri:) = o(/i"-*-i). (B.io) 

In other words, before ()B.4|) can be solved for h^^'' , it is necessary that the perturbing 
stress tensor be adjusted to be conserved with the metric g^^ and to satisfy ()B.10|1 . 
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Appendix C. V^T = implies the geodesic equation for a point mass 

We follow an example in reference jlH]. In 5^ (a;" — X°'{s))/ y/—g is a scalar field, 
and the factor u'' may be defined as a vector field by extension, in any smooth manner, 
away from the world line. Then, 



{g'a + U%a)VbT^' = fi{g\ + U%a) j 



oo 



(VfeM'')M' 



b 



-9 

b 



-9 

ds 



=^ r (Vbu^^n _ ^^^^^ 

J-oo 

where the second equality follows from properties of the projection operator g'^a + W^Ua- 
If VbT"'^ = 0, then it necessarily follows that the coefficient of the delta function must 
be zero for all proper times. A consequence is that v!'VhU"' = 0, the geodesic equation. 
A more formal proof of this result is in Poisson's review of the self-force ^S] , P 89. 

Appendix D. Gauge invariance of Eab{h) 

For a background geometry which is a vacuum solution of the Einstein equations, an 
infinitesimal gauge transformation, = x"" + with = 0(/i) changes the metric 
perturbation, = hab — 2V(a^b) + 0(/i^). But the operator Eabih) is invariant under 
such a coordinate transformation, 

^a6(VO=0. (D.l) 

This result follows immediately from the fact that the change in the perturbation of the 
Einstein tensor Eab under a gauge transformation is the Lie derivative of the background 
Einstein tensor £^Gab- For a vacuum background spacetime, this is zero. 
Equation ()D.1|) also follows from direct substitution into 

2Eab{h) = V^Kb + VaVbh - 2V(,V^/lf,)e + 'IRabKa + gab^^'^^K^ ' V'/i) (D.2) 

with hab = 2V(a^f,). It is easiest to consider the factor of gab separately, 
factor of gab = V^V^Ve^d + V^V^Vrf^c - 2V^VaV% 
= 2W''Vc^d - 2V'^Va 

= 0. (D.3) 

The second equality follows after use of the Ricci identity on the first two indices of the 
second term, use of Rab = for a vacuum spacetime and a relabeling of the indices. 
The final result follows after use of the Ricci identity on the second term of the second 
equality and use of Rab = for a vacuum spacetime. With hab = 2V{aC,b), the remainder 
of Eab{2V0 is 

remainder = Wc^a^b + ^'^c^bU + 2VaVbV% 

+ 2Ra''b''Vc^d + 2Ra%''Vdic. (D.4) 
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The analysis of this expression is lengthy but not difficult. It begins with using the Ricci 
identity upon the second and third indices of the ffist, second, fourth and sixth terms 
and upon the first and second indices of the fifth and seventh terms. The resulting 
terms with three derivatives may be paired up in a way to use the Ricci identity again 
and to reduce the entire expression to one involving only single derivatives. This also 
requires application of ()A.2|) . That the entire expression is zero, then follows from the 
symmetries of the Riemann tensor. 

Appendix E. Green's theorem for Eab 

The operator Eab{h) in (fTH|l . with an arbitrary tensor k°'^, satisfies the identity 

2k'''Eab{h) = V,F\k, h) - (F^ Kb). (E.l) 

where 

F%k, h) = k'^'V'^Kb - hv'h - 2(A;'=^ - ^g^k)V%hab - ^gabh) (E.2) 
and 

hab) = V'k^'Vchab - ^V'kVch 

- 2Va{k''' - ^g'''k)V\hbc - ^gbch) - 2k'''R,\''K,. (E.3) 

Note that the "inner product," {k"''', hab) = (^"^, kab) is symmetric under the interchange 
of and kab- It follows that 

k'^'Eabih) - K'Eabik) = [F%k, h) - k)] . (E.4) 

Which is a tensor version of Green's theorem for the differential operator Eab{h). 

The derivation of equation ()E.1|) is straightforward. Contract ()15p with an arbitrary 
symmetric tensor /c*^*, and move fc"^ inside Va in each term by "differentiating by parts." 
the divergence terms determine F'^{k, h). 

Appendix F. Singular gauge transformations 

Let be a, possibly distribution valued, vector field. And let hab = — 2V(a^f)), as for 
a gauge transformation. Also, let kab be a smooth "test" tensor with compact support. 
Then 

:.abrr> / u\ I / lab 



k'''Eab{h)4^gd^x= \ K'Eab{k)^gd^x 



2 JiV''e)Eabik)y^d\ (F.l) 



from ()E.4|) . after dropping the divergence term. An integration by parts and application 
of the perturbed Bianchi identity ()A.4|) yields 

J k-''Eab{h)^gd^x = 2 j [Eabik)]^d''x 

= 0. (F.2) 
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Thus, we demonstrate that given a solution to the inhomogeneous, perturbed Einstein 
equations (fT^ , even a distributional gauge transformation leaves a distributional valued 
metric perturbation that continues to satisfy the perturbed Einstein equations in this 
distributional sense. 

Appendix G. Black hole moving through an external background geometry 

When a small Schwarzschild black hole of mass m moves through a background 
spacetime, the hole's metric is perturbed by quadrupole tidal forces arising from 2-ffa6 
in ()24|) or ()53|). and the actual metric near the black hole, including the quadrupole 
perturbation, is 

= ^?af" + 2hat + o(rV7^^), (G.l) 

where the quadrupole metric perturbation 2hab is a solution of 

E!^-{,h) = 0. (G.2) 

Here -Ef^'^™ is the Schwarzschild geometry version of the operator given in (fT^ . The 
appropriate boundary conditions for ()G.2jl are that the perturbation be well behaved 
on the event horizon and that 2hab in the buffer region, where n <^r <^Ti. 

Our analyses of the boundary conditions and solutions of equation ()G.2|1 for slow 
motion are very strongly influenced by Poisson's recent analysis |27| |2H1 1211 of the same 
situation. In this appendix we describe 2hah up to a remainder of 0(r^/72.^). 

The appropriate boundary conditions at the future event horizon are that hah 
be finite and well behaved in a coordinate system which is well behaved itself. The 
Eddington-Finkelstein ingoing coordinates are satisfactory, and 

V = t + r^ and R = r, (G.3) 

where = r + 2mln(r/2m — 1); the angles 9 and remain unchanged. 

The odd and even parity parts of the metric perturbation are governed by the Regge- 
Wheeler jT^j and Zerilh [T7j equations, respectively. For our purposes, we change these 
equations from the Schwarzschild t,r to the Eddington-Finkelstein V, R. For example, 
the Regge- Wheeler equation for WiV, R) becomes 

where we have assumed that the angular dependence of W corresponds to a linear 
combination of £ = 2 spherical harmonics. We next assume that W is slowly changing 
in V and accordingly let 

W{V, R) = BWoiR) + B'Wi{R) + ... where B' = dB{V)/dV. (G.5) 

S is a function of V, and of the angles, that is related ultimately to the time-dependent 
external quadrupole moment of the geometry through which the black hole is moving. 
Thus B = 0{TZ~^) and B' = 0{TZ~^), in keeping with the requirement of slow time 
dependence. 
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With the form ()G.5|) substituted into ()G.4|) . we separate the terms of 0(7^~^) from 
those of 0{TZ~^) to obtain two equations. One is an ordinary, homogeneous differential 
equation for Wo{R), and the second is for for Wi{R) with a source term from the 
d^{BWo)/dVdR term. An analytic solution of the first equation is Wo{R) = R^. The 
solution of the second for Wi{R) is also analytic but more complicated, and constants of 
integration may be chosen so that Wi is well behaved at the horizon. This procedure thus 
provides a general solution for W^(V, R), up to a remainder of 0(7^^^). The even parity 
Zerilli equation may be solved in a similar manner or by using the simple relationship 
between solutions of the Regge- Wheeler and Zerilli equations |3U] . 

From the solutions of the Regge- Wheeler and Zerilli equations, the actual metric 
perturbations are determined by taking derivatives of the master variables. This results 
in a metric perturbation, in the Regge- Wheeler gauge, whose non-zero Schwarzschild 
components, as functions of V and r, are 

= - (r - 2m)2[f - 2mln(r/2m) S'^^'^] 

+ ^(3r^ - 12r^m + SGmV^ - 16mV - Sm^)^'^^) (G.6) 



^ r(2r^ — 3mr^ — 6m^r + 6m'^)^,(2) 

3(r — 2m) 



(G.7) 



l/i*- = - (r^ - 2m2)[^(2) _ 2mln(r/2m)r(2)] 

+ — (3r^ - ISmV^ - 12mV + Sm^)^'^^) (G.8) 



3r 



(2)1 



hrr = -r^lS'^^^ -2mln{r/2m)S' 

1 (3r^ — 12r^m + SGm^r"^ — 16m^r — 8m^) ,(2) /n 

+ 3 (7^2mj'^ ^ ^^-^^ 

= ^r{r - 2m)[i35f^ - 2mln(r/2m) Sf^] 

- ^(3r^ - 6r^m - 12rW + 12rW + 8rm^ + 8m^);B;J^^ (G.IO) 



9r 



,4 



This metric perturbation was first derived by Poisson [201 in a different gauge. 

In these expressions B^'^^ and ^^^^ are ^-dependent linear combinations of the i = 2 
spherical harmonic functions Y2^m{(^,4>)- The V and R coordinate components are all 
well behaved on the future event horizon. The V^-dependence of S^'^^ and B^"^^ shows 
that the metric perturbation propagates toward the black hole from a great distance as 
expected. 

To make contact with the actual, external geometry it is useful to expand the 
expressions given in (jG.Gjl - fjG.llj) for r in the buffer region, where m <^ r <^7l, and we 
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take advantage of the fact that B'^'^^ and S'^'^^ are 0(7^ ^). Thus, for TZ a Taylor 
series about V = t provides 

B{V) = Bit + r,) 

= S(^)+r,^ + 0(7^-^). (G.12) 

For m <^ r TZ, the even parity part of the metric perturbation in Schwarzschild 
coordinates is 

h2 dx" dx^ = - ^(2) [{r - 2mf df + dr^ + (r^ - 2rr?)aAB dx^ dx^] 

H \2{r + m) dt^ + 2(r + 5m) dr^ + (2r + 5m)(jAB dx^ dx^l 

_ /(2r^-3mr^-6mV + 6m3)^(,) ^^^^ ^ 0{m^£^^ /r% (G.13) 
3(r — 2m) 



and the odd parity part is 

'-{r-2m)Bf+^^{3r + Am)Bf 



hit dx" dx'' 



dt dx^ 



4 

+ 2—^ -ijj^ dr dx^ + 0(m^B^^^/r^). (G.14) 

12 (r — 2m) 

In this form E^"^^ and are considered functions of t and S^'^^ denotes the t derivative 

of 
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